Magnetic loops in the quiet Sun 
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ABSTRACT 

We investigate the fine structure of magnetic fields in the atmosphere of the quiet Sun. We 
use photospheric magnetic field measurements from SUNRiSE/IMaX with unprecedented spatial 
resolution to extrapolate the photospheric magnetic field into higher layers of the solar atmo- 
sphere with the help of potential and force-free extrapolation techniques. We find that most 
magnetic loops which reach into the chromosphere or higher have one foot point in relatively 
strong magnetic field regions in the photosphere. 91% of the magnetic energy in the mid chro- 
mosphere (at a height of 1 Mm) is in field lines, whose stronger foot point has a strength of more 
than 300 G, i.e. above the equipartition field strength with convection. The loops reaching into 
the chromosphere and corona are also found to be asymmetric in the sense that the weaker foot 
point has a strength B < 300 G and is located in the internetwork. Such loops are expected to 
be strongly dynamic and have short lifetimes, as dictated by the properties of the internetwork 
fields. 

Subject headings: Sun: magnetic topology — Sun: chromosphere — Sun: corona — Sun: photosphere 



1. Introduction 

The Sun's magnetic field lies at the heart of the 
heating of the solar corona and the solar chromo- 
sphere (with t he possible exception of th e basal 



flux; see, e.g., iBello Gonzalez et alJlioioh . Un 



fortunately, however, the magnetic field is mea- 
sured almost exclusively in the solar photosphere 
and needs to be extrapolated from there in or- 
der to obtain its structure in the Sun's upper 
atmo sphere. The balloon-borne Sunrise mis- 
sion (jSolanki et alj|20ld: iBarthol et al.ll2010l ) ob- 



tained the magnetic field in the quiet solar pho- 
tosphere with a very high and homogeneous spa- 
tial resolution. This allows us to investigate the 
3D structure of the quiet Sun's magnetic field 
in more detail as compared with SOHO ( e.g. 
Wieeelmann fc Solankil 12004 iHe et all l2007h or 



Hinode fe.g. | Martinez Gonzalez fc Bellot Rubio 



2009L iMartmez Gonzalez et all^oiof T 

The resolution of the extrapolated field in the 
vertical direction depends on the spatial scales of 
the photospheric measurements. With a pixel size 
of 40 km on the Sun and a spatial resolution of 
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~ 100 km, we can, for the first time, resolve the 
thin layer of the photosphere and the lower chro- 
mosphere with several grid points. This allows 
us to study the magnetic connectivity between 
photosphere, chromosphere and corona by well- 
resolved magnetic loops. We also briefly discuss 
the possible implications of our results for quiet 
Sun loop heating models. This is of particular 
interest due to the large number of horizontal 
magnetic features in the photosphere, many of 
them associated with emerging magnetic loops 



(Martinez Gonzal ez et al.l 12008c iDanilovic et al 



2010l ) and t he lack of acou s tic w ave energy flux 



proposed bv ICarlsson et al. ( 2007 ) which inspires 
us to co nsider alternatives, a lt hough the recent 
work by iBello Gonzalez et all (|2009l |2010| ) sug- 
gests that the acoustic energy flux may well have 
been underestimated in the past due to insufficient 
spatial resolution. 

2. Extrapolation of photospheric magnetic 
field measurements into the upper solar 
atmosphere 

Here we use a (phase diversity reconstructed) 
Stokes vector map of a 37 x 37 Mm quiet Sun 
region in t he photosphere recorded by Sun- 
RlSE/IMaX (jMartinez Pillet et al.ll2010l the data 
set was observed for 1.616 hours starting at 00:00 
UT on 2009 June 9th.). The magnetic vector is 
obtaine d by inverting these d ata using the VFISV 
code bv lBorrero et al.l ( 2010h . We take 200 km as 
the average height at which the Fe 1 spectral line 
at 5250.2 A senses B, although it is difficult to 
define ( Sanchez Almeida et al. 19961 ) and, in fact, 
varies from point to point over the solar surface. 
The height of this layer corresponds to z = in 
Fig.ffl ' 

To get a first impression of the 3D magnetic 
field structure in the chromosphere and corona we 
extrapolate the photospheric measurements into 
the atmosphere under the force-free assumption. 

VxB = aB (1) 
VB = (2) 
B = B bs in photosphere, (3) 

where B is the magnetic flux density and a is 
the force-free parameter. Here we use mainly po- 
tential field extrapolations (a — 0) due to the 



generally low signal level in the quiet Sun, in par- 
ticular of the linear polarization. Taking into ac- 
count the non-force-free character of photosphere 
and lower chromosphere (as, e.g., proposed in 
Petrie fe Neukirch 20001 IWiegelmann fe Neukirch l 



20061 ) require information on the plasma density 
and temperature distribution. The employed ap- 
proximation is su pported by spectro-polar i metri c 
observations from lMartmez Gonzalez et al.l ( 2010h , 
which suggest that quiet sun loops show a poten- 
tial field like structure. We solve equations ([T])- 
(03) with the help o f a fast Fourier approach (see 
Alissandrakid 198ll ) with the measured vertical 
magnetic field as boundary condition. A Fourier 
approach is justified because the magnetogram 
is almost flux balanced with ^2B Z /^2\B Z \ = 

x,y x,y 

—0.077. The computations are carried out in a 
cubic box with 936 x 936 x 468 grid points in 
x,y,z, where z is the vertical direction. We use 
a constant grid resolution of 40 km in each di- 
rection, which corresponds to a total 3D model 
volume of 37 x 37 x 19 Mm. Figure QJi shows 
the line-of-sight photospheric magnetic field and 
slices of the reconstructed magnetic field at differ- 
ent heights in the solar atmosphere in Fig. [Tp-d. 
With increasing height in the solar atmosphere 
the magnetic field becomes smoother and shows a 
dipolar structure at coronal heights. 

3. Statistics of magnetic loops. 

We compute magnetic field lines from all pixels 
above a certain threshold value, here \B Z \ > 30 
G in the photosphere, which corresponds to the 
3tr noise level of the magnetogram. The field line 
integration is started in a rectangular area at the 
center of the magnetogram, 150 pixels or 6 Mm 
away from the lateral boundaries. A total of 28005 
field lines are followed from one foot point to the 
other. Of these 1936 (7%) leave the computational 
domain through the lateral or top boundaries. We 
cannot say if these lines are really open or close 
outside the SuNRiSE/IMaX FOV. Consequently, 
we do not consider them further in this section. 
For simplicity we call closed magnetic field lines 
loops for the rest of this paper. Such magnetic 
loops are ubiquitous in the solar atmosphere, al- 
though only a fraction of these loops might be vis- 
ible in images made in chromospheric, transition 
region or coronal radiation. Loops with a height 
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Fig. 1. — Vertical magnetic field at different heights. Panel a) shows the measured magnctogram in the pho- 
tosphere and panels b-d) the vertical magnetic field distribution in different heights in the solar atmosphere, 
as reconstructed with a potential field model. The colour bars show the magnetic field strength in G. 
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Fig. 2. — a) A random selection of 2% of magnetic 
loops, computed from a potential field reconstruc- 
tion, b) Loop statistics for all resolved loops. The 
solid histogram-style line corresponds to a poten- 
tial field and the dotted and dashed lines to a lin- 
ear force-free model with a ■ L = 4 and — 4, respec- 
tively, c) Scatter plot of loop height vs. leading 
foot point vertical field. 



of less than 100 km (9321, 33% of all field lines), 
may not be spatially sufficiently resolved and are 
not considered further, leaving 16748 field lines 
(60% of all) to be studied. 

Figure [2^, shows a randomly chosen fraction of 
2% of these loops and Fig. [2J) a histogram of the 
loop height distribution for all resolved loops. The 
average loop height is 1.24±2.45 MmQ. The num- 
ber of small loops (below about 1 Mm) is signif- 
icantly larger than that of higher loops. 51% of 
the resolved loops d close within phototospheric 
heights (< 500 km) and 29% in the chromosphere 
(500 - 2500km). Finally 20% of the loops reach 
into the corona (> 2500 km), but half of these long 
coronal loops do not close within the Sunrise FOV. 

We also investigated how using a linear force- 
free model instead of a potential field model in- 
fluences the loop height statistics. These mod- 
els contain the force-free parameter a, which can- 
not be deduced from the available data. Conse- 
quently one unique linear force-free model cannot 
be computed. To investigate the influence of lin- 
ear electric currents, we carried out computations 
with aL = ±4, respectively, where L = 37 Mm 
is the size of the employed portion of the magne- 
togram (see dotted and dashed lines in Fig. [2p, 
respectively). 4 is already quite a large value for 
the normalized force-free parameter aL and close 
to the mathematical allowed maximum val ue of 
\aL\ = \/2k (see discussion in lSeehafer Il978l for a 
suitable rank of a- values for linear force-free mod- 
els). As visible in Fig. [2J) the loop height statis- 
tics are almost identical to those obtained from 
potential fields. Small differences occur only for 
field lines higher than about 5Mm. For our study, 
which mainly concentrates on lower chromospheric 
loops, it is therefore justified to consider only po- 
tential fields. 

Figure EJ; contains a scatter plot of loop height 
vs. magnetic field strengths at the leading foot 
point, i.e. the foot point with the largest field 
strength. It shows that loops of any height can 
begin from weak fields, but only long loops start 
from strong-field foot points. This means that 



x For higher (lower) threshold values of \B Z \ the number of 
small loops decreases (increases), which affects the average 
loop height, e.g., to 1.38Mm and 1.09Mm for a threshold 
of 40G and 20G, respectively. 

2 All loops higher than 100 km, including the ones not closing 
within the Sunrise/IMaX FOV. 



the stronger field (|£? 2 |>800G) loops are all long 
enough to reach into the solar corona, whereas the 
apex of a loop originating in a weaker field re- 
gion can lie in the photosphere, chromosphere or 
corona. 

4. On the origin of chromospheric mag- 
netic fields and magnetic energy 

In the following we investigate in more de- 
tail the origin of magnetic fields in the mid- 
chromosphere at a height of f Mm and are in 
particular interested in the magnetic connectiv- 
ity with photospheric fields. To this aim we start 
the magnetic field line integration in the chro- 
mosphere (see Fig. [TJd) at every pixel except in 
a layer of 150 pixels towards the magnetogram 
boundaries. We track all field lines down to the 
photosphere and determine here in particular the 
leading, stronger photospheric foot points E|. From 
this foot point map we can then distinguish be- 
tween photospheric pixels hosting magnetic loops 
reaching into the chromosphere or higher and pix- 
els not hosting such loops. Figure |3K shows the 
magnetic field strength distribution \B Z \ in the 
photosphere. The solid line is for all pixels, the 
dashed line for pixels hosting foot points of field 
lines that reach to a height of at least 1 Mm 
and the dotted line for pixels not hosting such 
loops. Above 500G the solid and dashed lines are 
practically identical and nearly all pixels above 
500G host footpoints of chromospheric loops. For 
weak fields (< 100G) only about 2% of the pixels 
host chromospheric loops, but the total number 
of these (hosting) pixels is still more than one 
order of magnitude higher than strong field pix- 
els. This raises the question, which photospheric 
fields can influence the chromospheric and coro- 
nal gas more strongly, the few regions with strong 
fields or the many regions with weak fields? To 
answer this question we compute the magnetic 
energy at a height of 1 Mm and investigate which 
fraction of the energy is magnetically connected 
to photospheric foot points with a particular field 
strength. As one can see in Fig. [3Jd the weak 



3 Some long field lines do not close within the Sunrise- FOV 
and for them we can only identify one foot point. For 
the analysis in this section these field lines are included in 
the statistics assuming that the known foot point is the 
stronger one. The relation to the weaker foot points is 
investigated in section [5] 
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Fig. 3. — a) Photospheric field strength \B Z \ 
statistics of the photospheric magnetogram (ex- 
cluding the 150 pixel layer towards the bound- 
aries). The solid histogram-style line corresponds 
to all photospheric pixels, the dashed line to pixels 
hosting foot points of field lines reaching at least 
to 1 Mm, i.e. into the mid-chromosphere, and the 
dotted line to pixels not hosting such field lines, 
b) Histogram showing the fraction of magnetic en- 
ergy at a height of 1 Mm (mid-chromosphere) in 
loops with a particular field strength of their lead- 
ing foot point. 
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photospheric fields (< 100G) hardly contribute to 
chromospheric magnetic energy. In addition, al- 
though only 0.91% of all photospheric pixels have 
a field strength above 100 G, they are connected 
to 97% of the chromospheric magnetic energy. 
The 0.32% of pixels above the equipartition field 
strength of 300 G contribute to 91% of the chro- 
mospheric magnetic energy. At the equipartition 
field strength (which is i n the range of 200-4 00G 



in the photosphere, see ISolanki et al.l 119961 ) the 



magnetic energy is identical to the energy of the 
convective flow. The equipartition field strength 
has relevance for quiet Sun magnetic features, be- 
cause magnetic flux tubes need to have a larger 
field strength in order to survive as a distinct 
feature. It is therefor natural to distinguish net- 
work elements from internetwork features by the 
equipartition field strength. Here we used the 
average value of 300G. 

Note that the photospheric field strength values 
are based on single component inversions assum- 
ing spatially resolved fields and are consequently 
lower limits (for the stronger field, th ey are close 



to the true value, see iLagg et al-lEoToh . and so are 



the above energy fractions of the chromospheric 
magnetic energy. The above numbers refer to the 
leading (stronger) foot point. 

5. Relation between the two foot points of 
magnetic loops 

In the last section we learned that 91% of 
the mid-chromospheric magnetic energy is topo- 
logically connected to photospheric foot points 
with superequipartition field strength. However, 
still unknown is if both foot points of long chro- 
mospheric and coronal loops have similar field 
strength or if it can be very different. If the latter 
is the case this would also mean that the multi- 
ple field lines starting within the area of a strong 
magnetic pixel will end in many, possibly widely 
distributed weak-field pixels 0. In Fig. Hi we 
show a scatter plot of weak foot point strength 
vs. strong foot point strength. 0. 95% of all loops 



4 In order to fulfill V • B = we get the relation B±Ai = 

B2A2 for the foot point areas A± and A2 of magnetic flux 

tubes containing a bundle of field lines. 
5 Necessarily, we had to exclude those field lines which do not 

close within the Sunrise-FOV, because we cannot identify 

both foot points for them. 
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Fig. 4. — Relation between strong and weak foot 
points of magnetic loops which reach at least into 
the chromosphere. Excluded have been loops that 
do not close within the Sunrise-FOV (14%), be- 
cause the strength of one foot point remains un- 
known. Panel a) shows a scatter plot of weak foot 
point strength vs. strong foot point strength. The 
solid diagonal line corresponds to equal strength 
of both foot points and the dashed horizontal and 
vertical lines mark the eqipartition field strength 
of 300 G. Panel b): Photospheric field strength 
\B Z \ statistics for the strong (solid line) and weak 
foot points (dashed line) of chromospheric and 
coronal loops. 
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with a leading foot strength of > 300 G, have the 
other foot point in regions with < 300 G. Figured)} 
shows the photospheric vertical field \B Z \ statistics 
for pixels hosting a leading (stronger) foot point 
(solid histogram-style line) and for pixel hosting 
a weak foot point (dashed line). Obviously the 
vertical field of both foot points for long loops 
(reaching at least into the chromosphere) is very 
different. Partly, this has to do with the differ- 
ent field strength distribution in the two magnetic 
polarities for this particular magnetogram. The 
maximum vertical field in the negative field region 
is 1311 G, compared to 707 G in positive polarity 
regions Q Investigations of further quiet Sun re- 
gions with high resolution are necessary to see if 
such different field strengths are common or an ex- 
ception (For quiet Sun investig ations of magnetic 



ception [t or quiet bun m vestig 
flux see, e.g. IWang et~ 19951 ). 



6. Conclusions and outlook 

We investigated the small-scale structure of 
magnetic fields in the solar atmosphere with em- 
phasis on the magnetic connectivity between the 
photosphere and chromosphere. We found that 
chromosphcric areas containing about 91% of the 
magnetic energy are topologically connected to 
a photospheric foot point with a strength above 
300 G, i.e. with an energy density higher than 
the average equipartition with that of connec- 
tive flows. For the majority of the loops the 
magnetic field strength in both foot points dif- 
fers significantly, with the second foot point hav- 
ing a field strength below equipartition, i.e. net- 
work elements connect magnetically mainly to in- 
ternetwork (IN) features, in g eneral agreement 
with n umerical experiments of ISchriiver &: Title 
(|2003l) ; IJendersie fc Peterl (|2006l) . An interesting 
question is to which extent the finding that net- 
work elements are connected mainly with IN fea- 
tures might i nfluence models of quiet-Sun loop 



heating (e.g Irlansteenl Il993l : IChae et al.l 12002 : 



Mulleret alll2003l . 1200 J ) . Such loop heating mod 
els usually assume a constant cross section, an ap- 
proximation which is not fulfilled for field lines 
connecting strong and weak field regions in the 



6 Indirect evidence of large network fields of 1-2 kG con- 
centrated in small r egions has been given already in 
iTarbell fc Title! l ll977f) . but the magnetograph resolution 
of that time has been too low to resolve such structures. 



photosphere. It is also worth noting that IN fields 
are more dynamic and short-lived than network 
fields (IN fiel d average lifetimes ar e about 10 min 



according to Ide Wiin et al.l 120081) . Since most 



chromospheric and coronal loops, at least in the 
observed region, are connected at one foot point 
with IN fields we expect the chromospheric and 
coronal field to be more dynamic than suggested 
by network fields alone. Time series of high reso- 
lution magnctograms in the quiet Sun can be used 
to investigate this further and to revisit coronal re- 
cycling (e.g., as investigated at much lower spatial 
resolution bv lClose et aDl2004l l2005h . 
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